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We show that the presence of primordial stochastic (hypercharge) magnetic fields before the elec- 
troweak (EW) phase transition induces isocurvature fluctuations (baryon number inhomogeneities) . 
Depending on the details of the magnetic field spectrum and on the particle physics parameters 
(such as the strength of the EW phase transition and electron Yukawa couplings) these fluctuations 
may survive until the Big Bang nucleosynthesis (BBN). Their lenghtscale may exceed the neutron 
diffusion length at that time, while their magnitude can be so large that sizable antimatter domains 
are present. This provides the possibility of a new type of initial conditions for non-homogeneous 
BBN or, from a more conservative point of view, stringent bounds on primordial magnetic fields. 



Large scale magnetic fields in diffuse astrophysical 
plasmas represent a well established observational fact 
since few decades. It has been realized through the years 
that magnetic fields coherent over diverse length scales 
at different epochs may have a variety of quite inter- 
esting phenomenological consequences. Magnetic fields 
coherent today over scales of the order of 30 Kpc are 
measured [jj] and have an important role in the dynam- 
ics of the galaxy, for example, in confining cosmic rays 
0. Magnetic fields at the nucleosynthesis epoch, even 
if not directly observable, could change the reaction and 
the expansion rate at that time. The success of the ho- 
mogeneous and isotropic BBN provides then interesting 
bounds on their existence || . 

There are neither compelling theoretical arguments nor 
motivated phenomenological constraints which could ex- 
clude the existence of magnetic fields prior to the nucle- 
osynthesis epoch. Moreover, to explain the origin of the 
galactic magnetic fields some authors often invoke the 
dynamo mechanism which might amplify the primordial 
"seed" magnetic field. It is a challenge to produce large 
scale seeds, and different ideas were aimed at this pur- 
pose. The energy scales involved vary from ~ 100 MeV 
for the QCD phase transition [J3J to ~ 100 GeV in the 
case of the EW physics MM and even closer to the Planck 
energy scale for inflation or string cosmology |7||^]. 

The purpose of this Letter is the study of the implica- 
tions of the primordial magnetic fields which existed even 
before the EW scale (i. e. for temperatures >100 GeV). 
The origin of these seeds is not essential for us and con- 
sequently we simply assume that they were generated by 
some mechanism before the EW phase transition. Our 
main point is that these fields produce baryon and lep- 
ton number inhomogeneities (isocurvature fluctuations), 
which could have an impact on the standard BBN. 

Let us start from some qualitative considerations. A 
unique property of "unbroken" U(l) gauge interaction 
is the absence of mass of its corresponding vector parti- 



cle. Static "magnetic" fields are never screened (in the 
absence of monopolcs) and thus homogeneous fields can 
survive in the plasma for infinite time. Under normal 
conditions (i.e. small temperatures and small densities of 
the different fermionic charges) the SU(2) xU(l)y sym- 
metry is "broken" down to XJ(1)em, the massless field 
corresponding to \J(1) em is the ordinary photon and the 
only long-lived field in the plasma is the ordinary mag- 
netic one. At sufficiently high temperatures T > T c , the 
SU(2)xU(l)y symmetry is "restored", and non-screened 
vector modes Y" M correspond to the U(l)y hypercharge 
group. Hence, if primordial fields existed at T > T c , 
they did correspond to hypercharge rather that V(1)em- 
There are essential differences between the interactions 
of magnetic fields and the ones of hyper-magnetic fields 
with matter. The ordinary electro-magnetic field has a 
vector-like coupling to the fermions, while the coupling 
of the hypercharge fields is chiral. Thus, if hyper-electric 
(£y) an d hyper-magnetic (7iy) fields are present simul- 
taneously, they cause a variation of the fermionic number 
according to the anomaly equation, d^j^ ~ j^jHy • £y 
(here g' the hypercharge gauge coupling constant). Now, 
the presence of non-homogeneous hyper-magnctic fields 
in the EW plasma with (hyper) conductivity er c always 
implies the existence of a related electric field, £y ~ 
—V x Tty Since for a general stochastic magnetic back- 

ground ((7xy • V x Hy) 2 ) ^ 0, the non-uniform hyper- 
magnetic field must produce baryon and lepton density 
perturbations because of the anomaly equation. In what 
follows we compute the amplitude of isocurvature fluc- 
tuations induced by this mechanism and discuss their 
physical relevance. 

The starting point of our discussion will be the general- 
ization of the magneto-hydrodynamics (MHD) equations 
( valid for ordinary electro-magnetic plasmas) to the case 
of hyper-magnetic fields with anomalous coupling to the 
fermionic degrees of freedom (see also (6j). These equa- 
tions have to be used for T > T r . We are interested in a 
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slow dynamics and we then assume that most of the parti- 
cle reactions are in thermal equilibrium in the expanding 
Universe (the list of those include all perturbative strong 
and weak processes, strong and EW sphalerons, Yukawa 
interactions of /i, r and s, c, b, t quarks). The parti- 
cle physics processes crucial for our purposes are those 
related to the U(l)y anomaly and to the slowest per- 
turbative reactions with right electron chirality flip (e.g. 
eji+Higgs — > e^+W). Thus, our variables are the space- 
dependent hyper-magnetic and electric fields Hy, £y 
and right electron chemical potential /ir^x). The gen- 
eralized Maxwell equations in a Friedmann-Robertson- 
Walker metric with scale factor a(r) are 

8H Y - - dE Y - -. 

— — = -V x E Y , -^r- + J Y = ■luaHy + V x H Y 

OT OT 

V-H Y =Q, V-E y =0, J y = ctE y , a(r)dT = dt (1) 

(Ey = <i 2 £y\ Hy = a H Y ; Jy — a3 jY] °~ — c c a). A new 
term, proportional the right electron chemical potential, 
comes from the anomaly contribution to the effective La- 
grangian of hypercharge gauge fields S , 



(2) 



Since the EW plasma conductivity is large, o~ c ~ o~qT 
with <to ~ 70 - 100 @, the time derivatives of the 
electric fields in Eq. ([!]) can be neglected (in the MHD 
context this is known as "resistive" approximation [flTf). 
This observations allows to express the induced electric 
field in terms of the magnetic one, 
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4 nH Y + V xHy 



(3) 



and derive an equation for Hy only. It is interesting to 
note that the presence of the fermionic chemical potential 
induces an electric field parallel to the magnetic one. 

The set of Eq. ([!]) has to be supplemented by the 
kinetic equation for the right electron chemical poten- 
tial, which accounts for anomalous and perturbative non- 
conservation of the right electron number: 



d_(VR\ = 9 12 783 
dt\T J 4n 2 cr c T 3 88 

-(r + r„)f , 



Hy-V xHy 



where T is the chirality changing rate, 
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(the numbers 783/88 and 783/22 come from the relation- 
ship between en number density and chemical potential 



S). An interesting consequence of Eqs. (^,||) is that in 
the presence of non-zero uniform magnetic field the right 
electron number is non-conserved (if T = 0), even for an 
abelian anomaly (cf. Ref. jL2|). 

Now we are ready to compute baryon number fluc- 
tuations produced in our scenario. We notice that at 
the temperature of the EW phase transition ~ 100 GeV, 
rt ^> 1. Then, since reactions with right electron chi- 
rality flip are in the thermal equilibrium, the adiabatic 
approximation can be used, and from Eq. (0) we have 
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Clearly, a non-uniform distribution of the right electron 
chemical potential induces baryon and lepton number 
perturbations of the same order of magnitude. We are 
not going to write the explicit formulae since there is an 
important "storage" effect which amplifies the estimates 
of Eq. (^) by many orders of magnitude. Equations 
pi) imply that 
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r + T H tr, 



■HyVxHy/0 
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Now, the change of Abelian Chern-Simons number is 
given by the time integral of (Q) . At the EW phase tran- 
sition the hyper-magnetic fields are converted into ordi- 
nary magnetic fields. The latter do not have coupling 
to the anomaly. Thus the CS number has to be trans- 
formed into fermions according to Eq. (Q). Inserting the 
coefficients, we arrive, from Eq. (R) at our main result: 
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(jib and s are the baryon and entropy densities, s = 
-g.TT 2 N e ffT 3 , N e ff is the effective number of massless 
degrees of freedom [106.75 for minimal standard model], 
M = Mpi/im^tff ^ 7.1 x 10 17 GeV). Notice that in 
Eq. (||) there is an enhancement by a factor ~ TMq/T 2 
arising from the time integration of the anomaly term. 
Some comments are now in order. 

(i) For the correctness of Eq. (||) the EW phase transi- 
tion should be strongly first order. Moreover, a necessary 
condition for EW baryogenesis [jl3| must be satisfied. In 
the opposite case all baryon number fluctuations will be 
erased by SU(2) sphalerons as it happens in the minimal 
standard model (MSM) fyj] , while this is not necessarily 
the case for the supersymmetric and other extensions of 
the standard model Q . 

(ii) Besides the primordial hyper-magnetic field, an es- 
sential quantity which fixes the amplitude of the isocurva- 
ture fluctuations is the rate of perturbative right electron 
chirality flip, T. For r# >r the amplitude of baryon num- 
ber fluctuations does not depend on the magnitude of the 
magnetic field fluctuations and it is determined just by 
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their spectral slope. For Th <r the rate of right electron 
chirality flip cancels out and the isocurvature fluctuations 
are fixed both by the magnitude and by the spectral slope 
of the primordial magnetic fields. In the MSM the rate 
r depends crucially upon the electron Yukawa coupling 
and is known to be T = Tt%, where T R ~ 80 TeV is 
the freezing temperature p6[ of the right electrons. This 
number appears to be too small to allow any interesting 
fluctuations. However, in the extensions of the standard 
model the rate T is naturally larger than in MSM. For 
example, in the MSSM the right-electrons Yukawa cou- 
pling is larger by a factor 1 / cos f3, which may increase 
the value of Tr by 3 orders of magnitude for experimen- 
tally allowed tan(/3) ~ 50. Cosmologically interesting 
fluctuations arise at T R > T* ~ 10 5 TeV @. 

We will assume now that r>r#, but similar conclu- 
sions hold true for the case T<r# and Tr > T* (for de- 
tails see ) . In order to compute the amplitude and the 
spectrum of the baryon number fluctuations we will also 
suppose that the Fourier modes of the magnetic fields are 
stochastically distributed, leading to a rotationally and 
parity invariant two-point function 



G i:j (r) = (Hii^Hjix + r)) , 



(9) 



where (...) denotes an ensemble average. In this case 
(S(^)(x,t))=0 Hi, but 

A(r, t c ) = J(8 (^) (x, t c )S (Pfj (x + f, t c )) . 

(10) 

Using the transversality of the magnetic fields it is use- 
ful to write the two-point function of Eq. (^|) in Fourier 
space 



G ij (k) = k 2 f(k)(S ij -^) 



(11) 



For /(fc) a power spectrum (modified by the typical ex- 
ponential decay of small scale magnetic fields given by 
Eq. (|l|)) is assumed 
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where k a 



k\ characterizes the strength of mag- 
netic fields, and e is the slope of the spectrum. A physi- 
cally realistic situation corresponds to the case in which 
the Green's functions of the magnetic hypercharge fields 
decay at large distance (i. e. e > in Eq. ([l2|)) and this 
would imply either "blue"( e > ) or "violet" (e ^> 1) 
energy spectra. The case of "red" spectra (e < 0) will 
then be left out of our discussion. The flat spectrum cor- 
responds to e < 1 and may appear quite naturally in 
string cosmological models Bj . 



The explicit result for Eq. (|1C|) at the EW phase tran- 
sition temperature is: 

a(r , tc) = J^^^ (1+0W) 
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(where T(z) is the Euler Gamma function and £ = 
k%/T c ). For a flat spectrum of magnetic fields (e <C 1) 
baryon number fluctuations may be rather large. For 
example, if the energy sitting in the background mag- 
netic field is comparable with the energy density of the 
photons, (H Y ) ~ ^ 4 t nen for the smallest possible scale 
r ~ l/fc CT ~ 10~ 9 x(EW horizon ~ 3cm) we get, from 

Eq. @ 8{n B /s) ~ j^jj^f ~ 10 3 . This number ex- 
ceeds considerably the measure of the baryon asymmetry 
of the universe ns/s ~ 10~ 10 , thus small size matter- 
antimatter domains are possible at the EW scale. At 
the same time, for even larger scales (possibly relevant 
for structure formation), the fluctuations of Eq. ( |l3| ) are 
quite minute (since their amplitude decreases with the 
distance as l/r 1+e ) and may be safely neglected. 

We consider now the question whether the fluctuations 
we found are able to affect the standard BBN. This de- 
pends upon the scale of fluctuations at T = T c . Short 
scale fluctuations (well inside the EW horizon) have dis- 
sipated by the nucleosynthesis time [[t9| through the 
combined action of neutrino inflation and neutron dif- 
fusion. Baryon number fluctuations affect BBN provided 
they are sizable enough over the neutron diffusion scale 
(3 x 10 5 cm) at the onset of nucleosynthesis (Tjys — 100 
Kev) (l9). The neutron diffusion scale, blue-shifted to 
T c ~ 100 GeV, becomes, L diff (T c ) = 0.3 cm. Taking 
again the flat spectrum for magnetic fields and assuming 
that their energy is ~ T 4 we obtain for the baryon num- 
ber fluctuations at that scale S(ns/s) ~ 10~ 5 ^> 10~ 10 . 
If magnetic fields are large enough, domains of mat- 
ter and antimatter may exist at the scales 5 orders of 
magnitude larger than the neutron diffusion length. Up 
to our best knowledge, there were no studies of non- 
homogeneous BBN with this type of initial conditions. It 
would be very interesting to see whether this may change 
BBN bounds on the baryon to photon ratio by changing 
the related predictions of the light element abundances. 
This possible analysis will not be attempted here. 

A more conservative attitude is to derive bounds on the 
magnetic fields from the requirement that homogeneous 
BBN is not spoiled, i.e. 



&(L diff ,t c ) < ^ . 



(14) 



In terms of £ and e (which completely define our stochas- 
tic magnetic background) the bound ( [H| ) becomes 
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FIG. 1. We compare the bound imposed by isotropy re- 
quirement at the BBN scale (Eq. ([Hi])) with the one imposed 
by homogeneity at the same epoch (Eq. @). We see that 
the homogeneity requirement (shaded region) is more con- 
straining then the one of isotropy. The plot is in the case 
Q B hl 00 = 0.01, h 100 = 0.6, a c /T c = 70. 



log£ < 



log §f -6.26+ \ log e + 14i 



log [fl B hl 00 ] 
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We plot it in Fig. [T| for a typical choice of the param- 
eters and 0.05<e<1.6 (for a c /T c ~ 70 - 100, fi B /i? 00 - 
0.1 — 0.01, 0.4 < /iioo < 1 this bound does not change 
significantly). The bound jl^ ) is quite strong for blue 
spectra (i. e. < e < 1). For violet spectra (e > 1 in 
Eq. (|l3|)) the fluctuations turns out to be parametrically 
smaller than (jib/s) at the neutron diffusion scale and 
then practically unconstrained by BBN. 

The bounds on ordinary magnetic fields at the nucle- 
osynthesis epoch do also apply in our case. For example 
in H was obtained that in order to be compatible with 
the isotropic nucleosynthesis Hn^ns) < 1 x 10 11 Gauss 
at a temperature T/vs ~ 10 9 °K which implies, at T c , 
-Hjv(£c) < 1.34 x 10 23 Gauss. By now comparing the 
magnetic hypercharge density with this bound we get a 
further condition in our exclusion plot, namely: 



log? < 



(11.30 - | log f?)e + loge- 0.2 



(16) 



This condition is reported in Fig. [I] (upper curve) and 
compared with the one of Eq. (|l5| ) (lower curve) . Accord- 
ing to Fi g. |l| Eq. ( |l6l) could be satisfied without satisfy- 
ing Eq. (|l5|) for e<1.4. This implies that the bound we 
derived is more constraining (by two orders of magnitude 
for magnetic field at e < 1) than the bounds reported in 
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Veneziano for helpful discussions and comments. 



[3. 

[4 
[5 

[6 
[T 

[8 

[9 
[10 



P. P. Kronberg, Rept. Prog. Phys. 57, 325 (1994). 
Y. B. Zeldovich, A. A. Rusmaikin and D. D. Sokoloff 
"Magnetic Fields in Astrophysics" (New York, Gordon 
and Breach, 1983). 

G. Greenstein, Nature 223, 938 (1969); J. J. Matese 
and R. F. O' Connel, Astrophys. J. 160, 451 (1970); D. 
Grasso and H. Rubinstein, Phys. Lett. B 77, (1996); B. 
Cheng et al., Phys. Rev. D 54, 4174 (1996); P. Kernan 
et al., Phys. Rev. D 54, 7207 (1996). 
G. Sigl, A.V. Olinto and K. Jedamzik, Phys.Rev. D 55, 
4582 (1997). 

T. Vachaspati, Phys. Lett. B 265, 258 (1991); K. En- 
qvist, P. Olesen, Phys. Lett. B 319, 178 (1993); T. W. 
Kibble and A. Vilenkin, Phys. Rev. D 52, 679 (1995); 
G. Baym, D. Bodeker and L. McLerran Phys.Rev. D 53, 
662 (1996). 

M. Joyce and M. Shaposhnikov, Phys. Rev. Lett. 79, 
1193 (1997). 

M.S. Turner and L.M. Widrow, Phys. Rev. D 37, 2743 
(1988); B. Ratra, Astrophys. J. Lett, 391, LI (1992); A. 
Dolgov and J. Silk, Phys. Rev. D 47, 3144 (1993); 
M. Gasperini, M. Giovannini and G. Veneziano, Phys. 
Rev. Lett. 75, 3796 (1995); Phys. Rev. D 52, 6651 
(1995); D. Lemoine and M. Lemoine, Phys. Rev. D 52, 
1955 (1995). 

A.N. Redlich and L.C.R. Wijewardhana, Phys. Rev. Lett. 
54, 970 (1985) 

M. Joyce, T. Prokopec and N. Turok, Phys. Rev. D 



53, 2930 (1 996); G. Baym and H. Heiselberg, astro 



ph/9704214 



[11] D. Biskamp "Nonlinear Magnetohydrodynamics" (Cam- 
bridge University Press, Cambridge 1994). 

[12] V. A. Kuzmin, V. A. Rubakov and M. E. Shaposhnikov, 
Phys. Lett. B 155, 36 (1985). 

[13] M. E. Shaposhnikov, JETP Lett. 44, 465 (1986); Nucl. 
Phys. B 287, 757 (1987). 

[14] K. Kajantie, M. Laine, K. Rummukainen and M. Sha- 
poshnikov, Nucl. Phys. B 495, 413 (1997) 

[15] M. Carena, M. Quiros, C. E. M. Wagner, Phys. Lett. B 
380, 81 (1996); M. Laine, Nucl. Phys. B 481, 43 (1996); 
J.M. Cline and K. Kainulainen, Nu cl. Phys. B 482, 73 
(1996); M. Losada, Eiep-ph/9605266 . 

[16] J. Cline, K. Kainulainen and K. Olive, Phys. Rev. Lett. 
71:2372, 1993; Phys. Rev. D 49, 6394 (1994). 

[17] M. Giovaninni and M. Shaposhnikov, in preparation. 

[18] If the spectrum is dominated by parity non-invariant con- 
figurations then, in general, (6 f^p) (S,t)) ^ 0. Exam- 
ples of this types of configurations can be found in Q. 
If this is the case, the baryon asymmetry of the universe 
can be a result of decay of Chern-Simons condensate Jl3| . 

[19] J. Applegate, C. J. Hogan and R. J. Scherrer, Phys. Rev. 
D 35, 1151 (1987); B. Banerjee and S. M. Chitre, Phys. 
Lett. B 258, 247 (1991);K. Jedamzik and G. M. Fuller, 
Astrophys. J. 423, 33 (1994); K. Jedamzik and G. M. 
Fuller, Astrophys. J. 452, 33 (1 995); H. Kurki-Su onio, 
K. Jedamzik and G. J. Mathews, a.stro-ph/9606011. 



4 



